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Transport of interacting electrons through a potential barrier: 
nonperturbative RG approach 
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PACS 71. 10. Pm - Fermions in reduced dimensions (anyons, composite fermions, Luttinger liquid, 
etc.) 

PACS 73 . 63 . Nm - Quantum wires (Electronic transport in nanoscale materials and structures) 
PACS 71 . 10 . -w - Theories and models of many-electron systems 

Abstract. - We calculate the linear response conductance of electrons in a Luttinger liquid with 
arbitrary interaction #2, and subject to a potential barrier of arbitrary strength, as a function of 
temperature. We first map the Hamiltonian in the basis of scattering states into an effective low 
energy Hamiltonian in current algebra form. Analyzing the perturbation theory in the fermionic 
representation the diagrams contributing to the renormalization group (RG) /3-function are iden- 
tified. A universal part of the /3-function is given by a ladder series and summed to all orders 
in #2- First non-universal corrections beyond the ladder series are discussed. The RG-equation 
for the temperature dependent conductance is solved analytically. Our result agrees with known 
limiting cases. 
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Electron transport in nanowires has been studied theo- 
retically for more than 20 years. Initially it was found that 
electron-electron interaction affects even the conductance 
of a clean wire [1,2] . In the case of realistic boundary 
conditions, namely attaching ideal leads to the interact- 
ing quantum wire, the two-point conductance of a clean 
wire is that of the leads, equal to one conductance quan- 
tum per channel, irrespective of the (forward scattering) 
interaction [3,4]. Alternatively, it has been argued that the 
screening of the external field by the interacting electron 
system leads to a renormalization of the conductance to its 
ideal value of unity [5] . The work of Kane and Fisher [2] 
and Furusaki and Nagaosa [6] showed, that interaction has 
a dramatic effect on the conductance in the presence of a 
potential barrier. For repulsive interaction these authors 
found that the conductance tends to zero as the temper- 
ature, T, or more generally, the excitation energy of the 
electrons approaches zero. This was shown at low tem- 
perature in the limits of weak potential barrier and strong 
potential barrier (tunneling limit), and at special values of 
the interaction parameter, K = \ and K = |, for all tem- 
peratures [2,7,8]. It has been argued, that these results 
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apply to a contact free four-point measurement, which is 
best realized by measuring the absorption losses of an a.c. 
field in the limit uo — > [3,4] (for a different view, see [5]). 

We recall that the reason for the strong suppression 
of the conductance by a repulsive interaction is that the 
Friedel oscillations of the charge density around the po- 
tential barrier act as a spatially increasingly extended ef- 
fective potential as the temperature is lowered. A proper 
treatment of the two-point conductance in the limit of 
weak interaction, taking into account the gradual build-up 
of the Friedel oscillations as the infrared cutoff is lowered 
has been given by Yue, Matveev and Glazman [9]. These 
authors used the perturbative RG for fermions to derive 
the conductance for an arbitrary (but short) potential bar- 
rier. A generalization of that approach to the case of two 
barriers has been given in [10]. 

In this letter we propose to extend the approach of Yue 
et al. to arbitrary strength of interaction. We argue that 
the /^-function of the RG-equation for the conductance 
can be obtained in very good approximation by summing 
a class of contributions in all orders of the interaction. As 
we show below, this is possible in this case, at least at low 
temperature, since the class of diagrams with the maxi- 
mum number of loops in any order, which are the ones 
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contributing to the /^-function, form a ladder series. At 
intermediate temperatures additional diagrams contribute 
small corrections to the /^-function. The result of the so- 
lution of the RG-equation for the conductance, using the 
approximate /^-function thus obtained, is found to agree 
with all known results on the scaling behavior of the con- 
ductance, including the case K = |, but goes far beyond: 
it is valid for any interaction strength K and any potential 
scattering strength. A careful analysis of the dependence 
of the result on the cutoff procedure chosen shows that cer- 
tain terms in the /^-function are not universal and depend 
on the cutoff scheme. 

The Model. We consider a one-dimensional system 
(coordinate x) of spinless electrons subject to a poten- 
tial barrier at x = 0. The barrier is characterized by 
transmission and reflection amplitudes t = t = cos#, 
r = — f* = isinOe 1 ^ with negligible energy dependence 
in the energy range of interest (width of order of temper- 
ature, T, around the Fermi energy, ep). We assume the 
extension of the barrier, a, to be narrow, akp <C 1, and 
neglect all interaction processes both close to the barrier 
\x\ < a (here k F is the Fermi wave number) and in the 
leads, \x\ > L. Beyond the large scale L the system is 
assumed non-interacting, which allows for an asymptotic 
single-particle scattering states representation. 

If c^ k and c 2k are operators creating electrons in right- 
moving and left-moving single particle scattering states of 
the barrier (k > 0), we may define the electron creation 
operator ^(x) as 

^(x) = {B(-x) [i/j+(x) + r^(-x) + #2>)] 

+ B(x) + r^(-x) + ^ 2 + (x)] } (1) 

where ^ 2 (x) = J °° ^e ±ikx c+ 2k . 

It is convenient to define current operators J M (x), \i = 
0,1,2,3 by J^(x) = ^ a ^ =1 ^i(ax)r^ (3 ilj ( 3(f3x), where 
r M are the Pauli matrices plus the unit matrix, -0+ (ax) = 
ipi(x) or i~ x ) f° r a = 1? 2. We call Jo the isocharge 
current and the vector J = (Ji, J2, J3) the isospin cur- 
rent; these operators obey U(l) and SU(2) Kac-Moody 
algebras, respectively [11,12]. In this paper we will not 
make use of these relations, but instead will work in the 
fermion representation. Nonetheless, the representation 
allows one to work in the chiral fermion representation. 

The particle density operators for incoming (i) and out- 
going (o) particles in terms of the J M 's are given by (here 
x > ) 



with R 33 = \t\ 2 - \r\ 2 = cos 20, R32 = Im{tP + tr*} = 
— sin 20 cos (/>, Rsi = Re{tr* — tr* } = sin 20 sin <p. 

We consider a model with interaction constant g<i (no 
backscattering, no Umklapp processes). The Hamiltonian 
is given by H = Hq + Hi , with 



Ho = 2ttvf I 



PiR,l(tx) = ^ 2 (=F^)^i,2(T^) = Jo(-x) ± J${-x) , 

Por(x) = [t^t{x)+f^t(-x)\ [t>i(a?) +f*il> 2 (-x)] 

= J (x) + Js(x) 

Pol(-x) = [ripl (x) +tip2(-x)] [r*ipi(x) + t*ip 2 (-x)] 

= J (x)-J 3 (x). (2) 



Here J3 = (RJ)s is the third component of the isospin 
current vector J rotated by the orthogonal matrix R^ v 



f 

Jo 



Hi 
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dx 



J 2 (-X)+J 2 (X)+J|(-X)+J3 2 (X) 

Jo(-x)J (x) - J 3 (-x)J 3 (x) , (3) 



where vp is the Fermi velocity. The forward scattering 
interaction of like-movers may be absorbed into redefi- 
nition of vp in the usual way [11]. In Fig. 1 we show in a 
pictorial way our parametrization of the fermionic densi- 
ties and the interaction, H\. In Fig. la the ^-interaction 
processes are shown in the usual scattering configuration. 
The representation in terms of the currents J M is in the 
chiral basis (all particles moving to the right, see Fig. lb), 
which leads to a seemingly nonlocal interaction. It should 
be clear that electrons in the left half space (x < 0) are 
not affected by the barrier yet, whereas electrons on the 
right (x > 0) are. Note that in the case of perfect reflec- 
tion, t = 0, we have J3 = —J3 and the observable densities 
(see below) form an Abelian U(l) sub-algebra of SU(2). 
This is the case of "open boundary bosonization" , which 
allows a complete and rather simple analysis. [13] One can 
also show that the J3 part of Eq. (3) can be reduced by a 
canonical transformation H' = HU to the Hamiltonian 
with the interaction part 

H[ = v F B.J(x = 0) - 2g 2 f dx J 3 (-x)J 3 (x) , (4) 

J a 

here U = expz f °° dx~B.J(x) and B = 2#(cos </>, sin </>, 0). 
[14] The first term in (4) corresponds to the rotation of the 
incoming isospin current J by the "magnetic field", B, at 
the origin, Fig. lb. Eq. (4) thus resembles the Hamiltonian 
for the Kondo problem in the current algebra approach. 
[12] The major simplification in our case is the classical 
nature of B, as opposed to the quantum Kondo spin S, 
see [12]. 

Current and conductance. The total electron density 
p(x) is given by 

p(x) = \pm(x) + Pol(x)} Q(-x) + [Por(x) + Pil{x)] Q(x) 
= Jo(-x) + J (a?) + sgn(x) [- J B (-\x\) + J 3 (\x\) 
= p c (x)+p s (x), (5) 



where subscript c(s) refers to isocharge (isospin) compo- 
nents. From the continuity equation d t p(x) = —d x j(x) — 
—i[p(x),H] one finds the current 



j(x) 



V F 



Jo(x) - Jo(-x) + J 3 (-\x\) + J 3 (\x\) 



jc{x) +js(x). 



(6) 



We now consider the linear response to an applied voltage 
V(x,i) = |y(t)sgn(x), which is seen to couple only to 
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Fig. 1: (a) The interaction of the left- and right-going densities 
in the basis of scattered waves, (b) The interaction of the left- 
and right-going densities in the chiral basis, corresponding to 
non-local interaction, Eq. (3). 



the isospin components, p s , of the density operator. The 
conductance is then given by (in units of e 2 /2irh) 



G(x,t) = -2nie(t) 



j 3 (x,t) , / dy p a (y,0) 
o 



(7) 



where we used the fact that correlation functions mixing 
the isocharge and isospin sectors vanish. 

Perturbation theory in g 2 . The contributions to G(Q m ) 
in n-th order of g 2 may be calculated with the help of 
Feynman diagrams in the position-energy representation 
(O m is the external Matsubara frequency). We draw n 
vertical wavy lines in parallel, the upper endpoint of the 
i-ih line at — X{ with isospin matrix , the lower one at 
Xi with matrix \R3ii T ap attached and carrying the factor 
—2g 2 ; ol, (/?) are isospin indices of ingoing, (outgoing) 
fermion lines. The external vertices are at — x with matrix 
|r3 and at y with matrix ^(Rr)s. The vertex points are 
connected by Green's functions 



Q{x ; u n ) 



v F 



-signK)9Kx)e-^^^ 



(8) 



where the uo n are Matsubara fequencies uo n = (2n + 1)ttT. 
All internal x -variables are integrated on the positive semi 
axis. The trace over the product of all isospin matrices in 
each fermion loop is taken and a factor of 1/n! is applied 
to each n-th order diagram. The limit Q m — >• +0 is taken 
at the end. 

The incoming component of p s (y), Js(—y), only con- 
tributes in zeroth order: G[ ^ = \. Adding the contri- 
bution from the outgoing component one finds 



G(°) = i(l 



cos 28) 



f 



(9) 



The diagrams of first order in 52 are shown in Fig. 2. 
Note that the "vertex correction type" diagram, Fig. 2c, 



(a) 





(c) 



Fig. 2: Three Feynman diagrams, depicting the first order in 
#2 contribution to the conductance. The other three diagrams 
are obtained from these by reversing the direction of fermionic 
propagation. 



is oc Q and will be dropped. One finds 



G (d 



— - sin (20) m — 
4tt v j T 



(10) 



in agreement with [9]. The ultraviolet cutoff, To, is de- 
termined by the width of the potential barrier a as To = 
i>ir/(27ra), the infrared cutoff arises at finite T > vpjL 
through the Green's function at time t = 0: Q(2a ; t = 
0) ~ Tj sinh(27rTa/^F)- In the limit of zero temperature 
the finite-T logarithm A = ln(To/T) is replaced by the 
zero temperature expression Ao = ln(T/a). 

In n-th order the diagrams with only one loop contribute 
the scale dependent terms [lnTo/T] n . Our prinicipal ob- 
servation here is that the diagrams with the maximum 
number of loops (n loops) contribute linearly in logarithm 
oc ln(To/T). They form a set of ladder diagrams. The di- 
agrams linear in ln(To/T) but not contained in this ladder 
series will be discussed below. The sum of all ladder dia- 
grams (see Fig. 3) may be calculated, and will be denoted 
by L(xi,X2]Wn). Later we will need the integrated quan- 
tity L(xi\uj n ) = J °° dx 2 e~ UJnX2 L(xi,x 2 ;uj n ), which 
obeys the integral equation 



L(x;u) = 



47r + ^(y + ^) / dze- uz L(z;uj) 



1 f°° 

+ -g 2 uj / dz e-^ x - z \L{z-uj). 

2 Jo 



(ii) 



Here we defined Y = cos 2$, g = ^2/27r and we use units 
with vp = 1. The Eq. (11) is of Wiener-Hopf type and its 
solution is L(x;u n ) = —A.ngbe~ UJnPX where p = yl — g 2 
and b = (1+p) / (1+p+gY) . The conductance contribution 
q( l ) i n linear order in In -p-, summed to all orders in g 2 , 
is given by 



G (L) 



l-Y 2 



dxidx2dyL(x±, x 2 \ uj)Q{x\ — x\ e) 



x G{-xi - x 2 ;e- co)Q(x 2 - y;e)g(y + X]e + Q) . (12) 



Taking the limit Q — > at x ^> vp /T one finds 



i + Vi 



g(l~^ 2 ) ln ?o 
2+gY T 



-5(1 - Y 2 ) In 



To 



(13) 



Comparing eq.(13) with eq.(10) we see that the resumma- 
tion corresponds to dressing of the interaction, g — » 
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X 2 



Fig. 3: Ladder series, L(xi, X2; cj n ), describing the combined 
effect of interaction and barrier, and leading to a linear-in- 
logarithm m(Xb/T) contribution to conductance. 

Renormalization group approach. In perturbation the- 
ory the n-th order contribution in #2 is a polynomial in A 
of degree n. If the theory is renormalizable, all terms with 
higher powers of A should be generated by a renormaliza- 
tion group equation for the scaled conductance, G(A), or 
equivalently for Y(A) = 2G(A) — 1. Our approximation to 
the /3-function of the RG equation is given by the prefac- 
tor of In t^t in the perturbation theory result (13), with Y 



replaced by Y(A): 

dY 
~dA 



2g{l -Y 2 ) 



1 + 



g 2 + gY 

Introducing the Luttinger parameter, K = [(1 
g)] 1 / 2 , we can rewrite (14) as 

dA 
dY ' 



Pl{Y). (14) 
*)/(! + 



1 



1 



(15) 



2(y + l)^" 1 - 1) 2(1 - Y)(l - K) 

which is easily integrated (see below). We note that (15) 
has a kind of duality symmetry: it is invariant under K — > 
if -1 , Y — > —Y. For weak interaction, we expand K ±x ~ 
1 =F g and recover the result by [9] . In the limiting cases of 
nearly transparent and nearly perfectly reflecting barrier, 
we recover the results by [2,6], with G ~ (T/Tb) 2 ^ -1 ^ 
and 1 - G ~ (T /T) 2(1_K) , respectively. 

A remaining question concerns the existence of addi- 
tional terms in the RG /3-function, not contained in the 
ladder series (15) (cf. [15]). To answer it, we have cal- 
culated G using computer algebra up to fourth order 
(~ 40,000 diagrams). [14] The results are summarized as 
follows. In higher orders we find both leading (~ g n A n ) 
and subleading (~ g n A k , k < n) contributions. Most 
of these terms correspond to Eq. (14). However, starting 
from third order of g we find also subleading contributions 
linear in A, which are explicitly different from the form 
(15) and arise from diagrams depicted in Fig. 4. Adding 
these terms to the /3-function of (14) we obtain 

dY 

— = -g(l-Y 2 )+g 3 (l-Y 2 ) 2 (c 3 -gYc 4 + 0(g 2 )), (16) 

with g the above result of the ladder resummation, C3 = 
7r 2 /12 and C4 = 0.238 — Notice the different power of 
(1 — Y 2 ) in front of the extra terms in (16), which renders 
these terms irrelevant in the limits G —> 0, G — > 1. 

Before proceeding further we should discuss the issue 
of universality of the logarithmic corrections, i.e. their de- 
pendence on the cutoff regularization scheme. Interpret- 
ing the ladder summation as a dressing of the interaction, 



t£3 c3 



+ symm. 



Fig. 4: Skeleton Feynman graphs, leading to lowest order log- 
arithmic contribution, g 3 A, beyond ladder series. 



Fig. 3, g — > g, the ladder result (14) thus corresponds 
to only one non-trivial Matsubara frequency summation, 
which amounts to a "one- loop" RG correction in the usual 
classification. In this case the scale invariant linear loga- 
rithmic contribution is not sensitive to the cutoff regu- 
larization, i.e. when L exceeds the inverse temperature 
InL/a — > ln(vi?/27rTa) and the prefactor of the logarithm 
is preserved. In other words, these contributions are uni- 
versal in the RG sense. In the order g 2 , non-trivial "two- 
loop" corrections to (14) are absent. In the order g 3 the 
"three- loop" contributions are divided into two groups. 
The first group consists of the first diagram in Fig. 4 and 
its symmetry-related partners: its contribution to G is 
7^Ag 3 (l — y 2 ) 2 ; this linear-in-A correction is again uni- 
versal. The additional diagrams in Fig. 4 contain both 
A 3 and A contributions. In this situation the linear-in-A 
terms contributing to the (3— function are dependent on 
the cutoff scheme, i.e. are non-universal. The above cited 
value C3 = 7r 2 /12 corresponds to T = , i.e. the hard in- 
frared cutoff L in real space. If we calculate corrections at 
T ^> vpL~ l , then we obtain the soft cutoff result C3 = 1/4 
instead, in agreement with [16] ; we use this latter value 
below. 

Inverting (16) we get dA/dY = (^(F))" 1 - g h(g,Y) 
where h(g,Y) = C3— gYO(\) + . . . for small g has been de- 
fined. Notice that in the non-ladder corrections g h(g, Y) 
we may substitute g by its renormalized value g — > g. 
Truncating at lowest level beyond the ladder series, we 
get 

dA/dY ~ [p^Y)]- 1 - c s g . (17) 



Integrating the latter equation we have 



(T/T ) 2(1_K) 
*(G) 



d>(G)/d>(G ), 
G K 



(18) 



1-G 



(K + G(l-K)) 



4c 3 (l-K) 



with C3 = 1/4 and the assumed initial condition G = Go 
at T = T . 

The equation (18) is the central result of this pa- 
per. Let us discuss it in more detail. The above ladder 
approximation(14) would correspond to setting C3 = in 
(18). It is seen that the above cited scaling-law depen- 
dences of G on T remain asymptotically exact at G — > 
and G — > 1. The existence of terms in the ^-function 



±u„ u „ — +. — u~ 



beyond the ladder series modifies the behavior of the con- 
ductance at intermediate values G ~ 1/2, and the role 
of h(g,Y) ~ 1/4 can hence be largely viewed as a re- 
definition of the cutoff energy In To — > In To + O(gGo), 
when going from higher to lower T. The duality symme- 
try, Eq. (15), between the scaling exponents is preserved, 
rp2(K-i) _^ rp2(K -i)^ - n con t ras t to the recent claim 
in [17] ; the breaking of duality reported in [17] might 
be connected to the approximate character of the solution 
for the set of flow equations there. 

Let us also compare our findings to exact expressions 
available from the thermodynamic Bethe Ansatz method. 
In the particular case of K = 1/2 {g = 3/5) the conduc- 
tance G is obtained as a closed analytic function of T. [2,7] 
Our Eq. (18) reduces to a quadratic equation, with rele- 
vant solution 

4r 2 T \t\^l-\\r\> 
1 + 4r 2 + VI + 16r 2 ' T T \ r \ 2 

(19) 

where \r\ are the bare transmission and reflection am- 
plitudes. 

At high temperatures the result reported in [2,7] tends 
to the clean limit value for a four-point measurement, 
G = |; here we assume that the two-point conductance 
is obtained by multiplying the latter result by a factor of 
1/K equal to 2 in order to make contact with our theory. 
Fixing the overall scale T* in both solutions (r = T/T*) 
by their high-temperature behavior, G ~ 1 — r _1 , we show 
the overall picture in Fig. 5. It is seen that the result of 
the ladder summation overestimates the renormalization 
of the barrier, predicting smaller conductance at low T. At 
the same time the adjustment of the ladder summation by 
C3 = 1/4, Eqs.(17), (18), (19), gives excellent agreement 
with the exact solution [2,7], with a relative deviation not 
exceeding 4 % in the whole temperature range. 

In case K = 1/3 (g = 4/5), relevant for the descrip- 
tion of point conductance between the quantum Hall edge 
states, [8] we do not have a closed analytic expression for 
G(T). It is however possible to make a comparison to our 
theory. We fix the overall temperature scale by adopting 
G ~ 1 — r -4 / 3 at high temperatures, then we have G ~ 9r 4 
from our Eq. (18). Multiplying again the result reported 
in [8] by a factor of 1/K = 3, we have G ~ 10.0638r 4 
at low T, whose pref actor is within 12 % from our value. 
We may thus conclude that our result (18) provides a very 
good approximation even in the strongly interacting case 
g ~ 1, i.e. it is sufficient for all practical purposes. 

Conclusion. In this letter we presented a theory of 
transport of interacting electrons through a potential bar- 
rier, in the linear response regime and at all temperatures, 
for any short range barrier and for any forward scattering 
interaction. We employed a representation in terms of 
chiral fermions, which greatly simplifies the perturbation 
theory in the interaction parameter g^. In this way the 
scale dependent contributions to the conductance may be 




T/T* 

Fig. 5: Conductance as a function of temperature for K — 1/2. 
The exact result available from TBA (multiplied by factor 2) 
is shown by pluses. [2,7] The results of ladder summation, (18) 
with C3 = 0, are shown by dashed line ; the result of Eq. (18) 
with the adjusting term C3 = 1/4 is shown by solid line. 

studied systematically. At low energies long range corre- 
lations lead to logarithmically diverging terms in the form 
of powers of A = In ^ . 

In particular, the terms linear in A may be summed 
to all orders in #2, and the prefactor may be identified 
with the /^-function of the renormalization group equation 
for the conductance as a function of the scaling variable 
A. At intermediate temperatures additional small correc- 
tions to the /3-function were found in third and higher 
orders of perturbation theory. Approximating these ad- 
ditional terms by the lowest order (in #2 ) gave excellent 
agreement with the known exact result at K = - . This ap- 
pears to be one of the few cases where the /^-function can 
be determined beyond perturbation theory. Our results 
are in agreement with all known results, where applicable, 
but go far beyond. The RG-equation may be integrated 
analytically to give the conductance as an implicit func- 
tion of the temperature. The method we describe here is 
quite general and may be of value for calculating transport 
properties of the model out of equilibrium or of other mod- 
els in which the /^-function may be obtained by summing 
a ladder series. 

* * * 
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